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1.  SUMMARY  OF  RESULTS  --  June  to  December  1972 


The  work  performed  in  the  period  of  June  1972  to  December  1972  under 
Contract  No.  F30602-72-C-0486  is  a  logical  extension  of  that  performed  under 
Contract  No.  F30602-71-C-0356  from  June  1971  to  June  1972.  The  technical  and 
final  reports  that  have  appeared  previously  under  both  these  efforts  are  listed 
under  [1]  in  the  References  and  are  cited  in  abbreviated  notation  (TRI  through 
TRVI)  in  the  text.  A  summary  of  the  work  reported  in  TRVI  and  of  subsequent 
work  follows. 

1.1  The  Furutsu  Theory 

We  have  indicated  that  Furutsu 's  result  for  the  statistics  of  the  log- 
irradiance  of  a  beam  wave  in  turbulent  air  may  be  incorrect.  The  error  appears 
to  lie  in  an  approximation  made  to  solve  the  equation  for  the  fourth-order 
coherence:  an  approximation  that  assumes  that  the  permittivity  correlation  may 
be  replaced  everywhere  by  the  square  law  (which  only  holds  for  very  small 
separation  distances) . 

1.2  The  XVII-th  General  Assembly  of  URSI 

We  attended  the  XVII-th  General  Assembly  of  the  International  Union  of 
Radio  Science  (URSI)  in  Warsaw,  Poland,  August  21-29,  1972.  A  special  report 
on  this  conference  has  been  issued  to  all  who  receive  Management  ReportB. 

Copies  are  available  upon  request. 

1.3  Resolution  Study 

A  simple  model  was  constructed  for  the  purpose  of  estimating  when  one 

can  no  longer  distinguish  two  bright  point  features  on  an  object  in  the  sky 

due  to  intervening  layers  of  turbulent  air.  The  criterion  of  clear  distinction 

between  two  such  objects  separated  by  a  distance  p  is  determined  by  the  ratio 

2  1/3  3  ° 

of  o  to  the  distance  C  K  '  L  .  When  this  ratio  becomes  small,  the  two  points 

n  m 

smear.  Details  are  worked  out  for  horizontal  propagation  and  for  slant-path 
propagation. 
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1.4  The  Log-Amplitude  Variance  (Plane  Waves) 

The  plane-wave  theory  of  TRV  has  been  completed,  and  a  calculation  of  the 
log-amplitude  variance  yields  <6X2>  -  O-dK^Vc2)'1'  In  the  saturation 
regime.  Comparison  with  plane-wave  data  appears  to  be  reasonable. 

1.5  Angle-of-Arrival  Power  Spectrum  for  an  Interferometer 

The  calculation  of  the  angle-of-arrival  power  spectrum  given  in  TRIV  has 
been  extended  for  the  case  of  a  simple  interferometer  of  two  parallel  rays. 

Some  closed-form  results  are  obtained,  and  -  specifically  -  the  single-ray 
formulas  are  shown  to  require  only  slight  modification  in  most  cases  of  inter¬ 
est.  An  important  difference  with  results  of  other  workers  appears  to  be  the 
prediction  that  W(u>)  -  uf2/3  in  the  inertial  subrange  of  turbulence,  and  not  > 

W(co)  “  uT8^3! 

1.6  Spherical-Wave  Irradiance  in  Turbulent  Air 

An  extension  of  the  theory  for  plane  waves,  mentioned  above  in  subsection 
1.4,  has  been  worked  out  for  spherical  waves.  It  yields  <«X  >  “  1*25(,cm  L  V 
in  the  saturation  regime.  It  is  the  topic  of  Section  2  of  this  report. 

This  completes,  and  even  extends  somewhat,  the  required  tasks  of  the 
effort  under  Contract  No.  F30602-72-C-0486.  Extensions  of  this  work  vili  be 
performed  under  a  future  contract  affott. 
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2.  IRRADIANCE  SCINTILLATION  FOR  SPHERICAL  WAVES 


The  purpose  of  this  section  is  to  expand  the  plane-wave  theory  of 
irradiance  scintillation  to  include  spherical  waves,  and  specifically  to  com¬ 
pute  the  variance  of  the  log-amplitude  in  that  case.  There  are  several  versions 
of  the  plane-wave  theory  -  in  diverse  stages  of  development  -  in  circulation, 
and  the  reader  may  have  difficulty  in  following  the  development.  A  brief 
guide  through  previous  work  will  therefore  be  given. 

The  first  of  our  reports  to  deal  with  the  subject  is  TRIII  (the  notation 
for  referring  to  previous  reports  is  explained  on  p.l).  It  Is  chiefly  a 
tutorial  work  that  develops  the  basic  equation  (6.10),  and  then  works  out  all 
of  the  theory  except  that  for  the  saturation  regime.  Sections  10  through  12 
of  TRIII  are  self-contained,  and  serve  to  prove  a  crucial  theorem  in  optical 
propagation  through  turbulent  air:  the  n-point  correlation  of  refractive- 
index  fluctuations  always  occurs  in  certain  integrals  such  that  it  may  be 
factored  into  one  ordered  product  of  two-point  correlations,  even  though  it 
is  not  a  multivariate  Gaussian. 

In  TRV,  a  first  attempt  was  made  at  a  saturation-regime  theory.  As  we 
see  it  now,  there  are  objections  to  be  raised  against  the  development  given 
In  Section  III  of  TRV.  Furthermore,  the  important  parts  of  the  theory  are 
really  in  the  Appendix,  which  has  been  very  abbreviated  and  which  yields 
Eq.  (A9)  [unfortunately  not  labeled  In  TRV]  as  s  result  of  a  development  from 
Eq.(AA)  on,  that  differs  somewhat  from  that  given  in  this  report.  Consequently, 
Sections  III  through  V  and  the  Appendix  of  TRV  have  been  thoroughly  revised 
for  publication  in  the  Journal  of  the  Optical  Society  of  America  in  February 

or  March  1973. 

Therefore,  there  appears  to  be  a  departure  from  a  fully  logical  develop¬ 
ment  in  Section  3  of  TRVI,  because  we  do  indeed  start  from  Eq.(A4)  of  TRV 
bearing  in  mind  that  the  development  of  Eq. (A4)  from  Eq.(A2)  has  been  greatly 
expanded  in  the  article  to  appear  in  the  J.  Opt.  Soc.  Am.  -  and  then  ignore 
the  further  development  in  TRV.  Section  3  of  TRVI  does  repeat  the  pertinent 
development,  leading  to  an  expression  for  the  log-amplitude  variance,  also 
gi/en  in  the  aforementioned  article. 
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Let  us  now  continue  the  development  for  spherical  waves.  In  referring 

to  TRIII ,  we  find  that  the  first  question  is  whether  a  fundamental  equation 

such  as  (6.10)  can  be  developed  for  spherical  waves,  using  the  Green's 

function  G  (r,r. )  of  Eq.(2.7).  When  this  Green's  function  is  used  in  Eq.(3.3), 
s  1 

we  can  define  a  stationary-phase  exponential  expCi^)  as  in  Eq.(6.1),  but  in 
this  case, 


V  -  k[(ApZ  +  Az2) 1/2  -  <p*  +  z*  1)A/^  +  (pA  h  zA)i,A]  -  K  *p  ,  (1) 

m  m  m  m-1  m-1  mm  mm 


1/2 


2x1/2, 


and  the  points  of  stationary  phase  are  found  by  setting  the  gradient  of  ¥ 

with  respect  to  p  equal  to  zero.  It  is  easily  seen  that  the  stationary- 
in 


phase  points  Pm  obey 


where 


Ap  /Ar 
m  m 


/k  +  p  /  r  , 
n  n 


(2) 


j=m 

We  adhere  to  the  notation  developed  in  previous  work:  it  is  quite  uniform 

and  consistent.  An  exhaustive  glossary  of  symbols  is  given  in  TRIII.  Note 

that  the  physical  content  of  Eq.(2)  is  obvious  and  practically  identical  to 

that  for  plane  waves.  The  initial  propagation  angle  is  not  zero  but  Pn/rn 

for  spherical  waves  because  the  source  is  a  point  at  p  B  0,  z  “  0,  and  the 

first  deflection  occurs  at  (p  ,z  ).  The  scattering  angli  at  r  is  given  by 
_►  n  n 

-Q  /k,  exactly  as  in  che  case  of  plane  waves.  Consequently,  the  physical 
m 

small-angle  approximations  are  identical  too. 

Since  this  identity  permits  one  to  utilize  Apm  «  Az^,  pffi_^  «  zm_^» 

and  p  «  z  in  Eq.(l),  one  can  approximate  V  adequately  by 
mm  01 


¥  » 
m 


k 

2 


z  z  , 
m  m-1 


(3) 


It  requires  some  algebraic  rearranging  of  terms  to  obtain  this  form,  but 

2 

once  we  have  Eq.(3)  we  can  carry  out  all  the  d  Pm  integrations  in  Eq.(3.3)  of 
TRIII  to  obtain 


Bn  -A  i azmJa  \  ‘[a\r  '  J 

n 

where  Q"  *  Sz  K 
m  j-m  J  J 

This  form  is  analogous  to  Eqs.(6.10)  and  Eq.(9.1)  of  TRIII.  As  in  the 
transition  to  (9.3),  we  prefer  to  rearrange  the  exponential  terms.  The 
ultimate  result  is 

Bn SI(YV  • 

'if*  '  r  Kmlexp  [- 7(L-zm)Bn,'^fl]  ’ 


where  Q’  ■  V^K.z  ,/L 
111  j"m  J  J 

This  form  is  the  spherical-wave  analog  of  Eq.(14a)  of  TRV,  which  is  the 

point  of  departure  for  nearly  all  results  (except  the  saturation-  egime) . 

Note  the  occurrence  of  the  extra  geometrical  factors  z  /L  which  make  these 

m 

exponentials  slightly  less  significant  than  in  the  plane-wave  case.  All  known 
results  for  spherical  waves  follow  from  Eq.(5).  For  example,  the  extremely 
important  intermediate  steps  discussed  in  Section  13  of  TRIII  yield  the  same 
result  as  Eq.(13.5)  for  the  bead  contributions,  i.e.,  for  correlations  between 
two  6e(K,z)  factors  of  one  and  the  same  Born  term  B^.  However,  correlations 

between  an  i-th  B  factor  and  a  j-th  B*  factor  of  <B(1) - B(N)B*(1) - B*(N)>, 

or  between  two  unasterisked  or  two  asterisked  factors  of  this  2N-point  corre¬ 
lation  ultimately  yield  Eq.(l6.8)  with  new  definitions, 

C  =  cos [K^(z  /L) (L-z  ) /k] 
in  mm  m 

cyij)  =  cob  fi-W^-Q^UL^/k] 


Note  that  the  Q'^  are  defined  in  Eq.(5)  above.  In  fact,  with  these 
definitions,  the  entire  development  of  TRIII  from  Section  16  on  can  be  con¬ 
tinued.  Likewise,  Eq.(18a)  of  TRV  also  can  be  taken  over  unchanged,  provided 

C i  and  C^(ij)  are  given  by  the  above  Eq.(6). 

The  modified-Rytov  result  for  spherical  waves  follows  as  before  by  setting 
c  (ij)  =  Thc  error  estimate  is  identical  to  that  developed  in  Eqs. (17.12) 
and  (17.13)  of  TRIII.  The  result  for  <6x2>»  the  log-amplitude  variance,  is 
identical  to  Tatarski’s[2]  Eqs. (49.25)  and  (49.36).  Equation (49. 36)  in  our 
notation  is 

<6X2>  =  0.124cV/6LU/6,  <7> 

different  only  from  Tatarski's  expression  because  C2  =  4C2.  The  Rayleigh 
result  is  unchanged  from  the  plane-wave  case  because  the  only  important  con¬ 
tributions  require  Q^(:i)«  Q^(i)  to  within  errors  of  order  l2/L  in  deviations 

of  C  (ij)  from  unity. 

However,  the  most  important  matter  is  to  derive  a  saturation-regime 
corollary  to  Eq.(7).  One  can  start  -  once  again  -  from  the  parabolic  equation 
for  the  normalized  field  B  =  rE  exp(-ikr): 


AtB  +  2ik(r £)B  +  k26eB  =  0,  (8) 

where  r  is  a  unit  vector  of  the  radial  coordinate,  which  is  not  going  to  de¬ 
viate  strongly  from  the  z-direction.  It  yields  for  =  InB: 

A  ijj  +  2ik  r*^  +  (k25e  -  2ik/r)  +  $T*P)2  =  0,  (9) 

which  yields  the  spherical-wave  Rytov  approximation  ^r(L)  when  the  last  non¬ 
linear  term  is  dropped.  I.e., 

li 

VL)  /-  exPti^2^L^_1  z(L-z)AT]5e(o,z).  (10) 

Returning  to  Eq.(9),  we  note  that  the  plane-wave  procedure  for  eliminating 
($  i>)2  can  be  repeated  verbatim.  This  procedure  was  not  given  explicitly  in 
TRV,  and  it  will  therefore  be  included  here,  based  on  the  J.  Opt.  Soc.  Am. 
manuscript. 
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We  integrate  Eq. (9)  along  geometrical-optical  rays  with  coordinates 
r  =  [p'(s),s]  where  s  is  a  ray  coordinate  and  ^pj  (s)  =  [x ' (s) ,y’ (s) ]  are  two 
mutually  orthogonal  coordinates  such  that  pi  (s)  and  s  are  the  three  components 
of  a  locally  orthogonal  coordinate  system.  Let  B(z)  be  the  angle  between  s 
and  z  after  traveling  a  distance  ^z  along  s.  Let  9^  and  9y  be  the  projections 
of  9(z)  upon  the  y  =  0  and  x  =  0  planes,  respectively.  It  can  be  seen  that 

9/9x  =  (9y/9)9/9x'  +  (9x/9) (l+92)~1/29/9y’  +  Bx(l+92)_1/29/3s 
9/9y  =  -(9  /9)9/9x'  +  (9  /9) (l+92)-1^29/9y'  +  9  (1+92)  L^29/9s 


9/9z  = 


-  (l+92)“1/29/9y’  +  (l+92)'1/29/9£ 


In  particular,  we  note  from  Eq.(ll)  that 
9/9z  -  9/9s  +  0(9) 


Ax  =  Ax*  +  0(9). 

Furthermore,  it  can  also  be  seen  from  the  coordinate  transformation  that 

9/9x'  =  (9  / 9)  9/9x  -  (9  /9) 9/9y  <13 

y  x 

9/9y'  =  (l+92)_1/2[(9  /9)9/9x  +  (9/9) 9/9y  -  99/9z]. 

x  y 

2 

With  Eqs. (11-13),  we  can  rewrite  Eq.(9)  after  dividing  by  k  , 

[2ik_19/9s  +  0(0) +  [k_2AT'  +  0(9) 

-2ik/r  +  6e(p'(s),s)  +  [(k  1V  '  +  0(9))^]2  =  0  . 

Now  we  shall  choose  9  and  thus  define  the  ray  coordinates.  To  do  so,  we 
define  the  complex-angle  vector 

l  =  -  ik'^X  =  -ik_iVT^  .  (1] 

Thus.  9  has  two  complex  components,  9  and  0  ,  defined  by  replacing^  in 
Eq.  (15)  by  x9/9x  and  y9/9y,  respectively.  The  magnitude  of  6x(0y)  is  the 
square  root  of  the  sum  of  the  squares  of  both  terms  of  the  9x(8y)  component 
of  Eq.(15).  Note  in  particular  that  92  =  k  2[(<Sq’(t))  + 
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Using  Eq.  (13)  with  Eq.(15),  we  see  that  =  0  witl  this  choice.  The 

last  term  of  Eq.(14)  is  therefore  zero,  except  for  corrections  of  0(0)  in¬ 
cluded  in  the  other  terms.  Consequently,  the  solution  of  the  thus  linearized 

Eq. (14)  becomes 


L 

w  =  f-jdz  exp  [i  (2kL)_1z(L-z) A^]  6e[p(z)  ,z] , 


where  we  have  replaced  ds  by  dz  and  L-s  by  L  z,  and  have  written 
6e(s)  =  6e(p(z),z]  in  terms  of  the  old  coordinates.  The  crucial  difference 
from  Eq . (10)  is  the  dependence  of  6e  upon  the  coordinate  p(z)  which  is  given 
by  the  fundamental  relationship,  for  z  £  L, 

dp/dz  =  0(z)  ,  (17 

with  $(z)  given  by  Eq.(15).  The  set  of  three  equations  (15-17)  comprises  a 
nonlinear  set  of  equations  for  ^(L).  The  procedure  for  solving  them  is  not 
quite  as  circular  as  it  may  seem.  Assume  that  we  know  \p(z)  for  z  <  L. 

Equations  (15)  and  (17)  then  determine  p(z)  for  z  <  L.  The  result  is  inserted 
into  Eq .  (16)  to  yield  ijj(L) .  At  the  very  least,  this  suggests  a  numerical 
procedure.  It  does  more  than  that,  however,  because  the  role  of  p(z)  appears 
only  to  be  a  description  of  the  refractive  properties  of  the  radiation,  not 
of  the  diffractive  properties.  It  is  for  that  reason  that  one  may  ignore 
diffraction  in  Eq.(15)  and  utilize  simple  ray  formulas.  Thus,  Eq.(15)  is 
approximated  well  by 


z 

8(z>  -  t0  +  \  fiz  i  ^66(0,5^) 

0 

z 

=  $  -  -i~  /dz.  /d2K  KSe(K,z.)  , 

o  q„2  J  1 


where  t  is  the  initial  angle  of  the  unscattered  ray  from  source  to  point  of 
scattering  with  the  z-axis,  and  6e(K,z)  is  the  two-dimensional  transverse 
Fourier  transform  of  as  usual. 

We  note  that  <Sx2(L)>=  -<X(L)>  where  the  mean  <x(L)>  is  determined  by 
the  real  part  of  Eq.(16).  It  can  be  seen  that  Eq.(16)  and  the  above  con¬ 
siderations  yield 


<ox  CL)> 


The  remaining  problem  is  to  substitute  Eq.(18)  into  Eq. (19)  and  perform 

the  calculation.  The  procedure  is  quite  comparable  to  the  plane-wave  case. 

The  first  step  is  to  integrate  from  z^  =  0  to  z^  =>  z  in  Eq.(18)  to 

obtain  p(z)  according  to  Eq.(17).  If  this  is  done  with  the  boundary  conditions 

p(0)  =  p(L)  =  0  one  obtains 
L 

P(z)  =  /dZl  M(L, z ,z^) J'd2KK&e(K,zi) 

0  (20) 

M(L,z,z^)  =  z^(L-z)/L  ,  Zy  <  z 

=  z(L-z1)/L  ,  z1>  z 

Two  steps  are  required  to  reduce  Eq.(19)  to  a  form  suitable  for  computing 
closed-form  expressions  [see  the  Appendix].  They  are 


<6e(K,z)  exp[-iK*p(z) ]>  =  <-iSe(K,z)K*p(z)>exp[-K2<p2(z)>/4] 
<-ifie(K,z)K*p(z)>  =  («2/4)K2$(K)  (L-z)z/L 


(21a) 

(21b) 


A  calculation  of  <p  (z)>  remains  to  be  performed.  We  utilize  Eq.(20), 
and  the  usual  approximations  as  in  the  Appendix  calculation  of  Eq.(A.8)  to  ob¬ 
tain  from  straightforward  calculation 

00 

<p2(z)>  =  («2/24tt)  /  dKK3$(K)  z2(L-z)2/L 

0  (22) 
=  4yr(7/6)C2<1/3  z2(L-z)2/l, 

where  y  is  defined  below  in  connection  with  Eq.(23).  The  remaining  work  is 
completely  analogous  to  the  derivation  of  Eq. (18b)  in  TRVI.  Equations (21a) 
and  (21b)  are  substituted  into  Eq.(19),  and  the  steps  described  beyond  Eq  ,  (18a) 
in  TRVI  are  repeated  verbatim  except  that  the  above  spherical-wave  forms  are 
used.  The  result  is 


/ 


1  00 

<5X2>  =  |yk7/6L11/6cjJdy  y5/6(l-y)2y*dx  x1/6  exp[-x02f  (y)  ] 

0  o 

where  y  -  0.033*2  *  0.326  end  f(y>  =  yr(7/6>yCl-y) 2-  After  performing  the 
dx  integration,  and  some  reordering  of  factors  we  obtain 


<«x2>  -  e^'Vc2)-1'6 


J. 

where  3  «  [Yr(7/6)  ]“1/6  *  \fdy  yrl/3(l-y)1/3 , 

0 

and  after  a  numerical  integration  we  obtain 


V>  -  i.zsc^W)-1'6  (25) 

Equation(24)  is  the  final  result  for  spherical  waves.  It  should  be  noted 
that  Eq. (23)  yields  the  result  <6y2>  «  0.005  L  Cn  in  the  limit  ^  ^  °»  a 
result  that  agrees  with  Tatarski's  [ref.  2,  page  251]  for  spherical  waves. 
There  is  therefore  as  much  measure  of  confidence  in  Eq.(25)  as  in  the  plane 

wave  result. 

Let  a2  «  0.124  k7/6L11/6C2,  the  well-known  spherical-wave  result  for  the 
log-amplitude  variance  in  the  Rytov  approximation  [see  Eq.(6)].  It  is  easily 
seen  from  Eq.(25)  that 

<6y2>  w  0.88(K2Lik)  7^36(<J2)  1^6, 

and  therefore  at  k  =  107  m  1  and  =  1000  that 


<6X2>  W0.  36 (a2)'1/6  at  L  =  1  km 


W0.41(a2)"1/6  at  L  =  0.5  km 

We  have  utilized  Eq.  (27)  as  a  basis  for  comparing  the  available  published 
data.  In  Fig.  1  we  have  plotted  data  by  Kerr[4],  Ochs [5],  and  Gracheva  et  al. 
[6].  Two  lines  in  the  saturation  regime  represent  the  mean  data  of  these  two 
sources.  The  standard  deviation  around  the  mean  is  of  the  order  of  50%;  thus, 
the  theoretical  curve  [Eq . (27) ]  appears  to  be  low  but  still  within  the 

standard  deviation. 
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Figure  1.  Log-amplitude  variance  <<$X^>  vs.  parameter  =  0.124  k?/6T 11/6^2 
for  spherical  waves:  data  and  theory. 

We  have  also  reexamined  plane-wave  data  by  Gracheva  et  al.  [6],  and  by 

Mevers  et  al. [7],  replotted  in  Fig.  2.  Data  spreads  are  of  the  order  of  50% 

around  the  mean.  The  theoretical  curve  is  based  on  Eq. (21)  of  TRVI,  and  it 

appears  to  be  systematically  on  the  low  side  (spherical -wave  data  appear  to 

agree  better  with  the  analytical  prediction) . 

The  spherical-wave  curve  appears  to  fit  the  data  slightly  better  than 

the  plane-wave  curve.  It  is  interesting  to  note,  particularly  from  Fig.  1, 

that  the  data  do  not  exhibit  an  outspoken  frequency  dependence  -  in  accord 

with  the  saturation-regime  prediction.  It  should  also  be  noted  that  the 

—7 /18 

theoretical  curve  is  proportional  to  <m  .  We  have  arbitrarily  chosen 

£  =  6  mm,  but  a  choice  of  £  =12  mm  will  move  the  curves  up  by  a  factor  1.3. 
o  o 
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APPENDIX 


Derivation  of  Equations (21) 

The  derivation  of  Eqs.  (21)  from  Eq.(19)  is  not  quite  trivial.  The 
correlation  <6e(K,z')exp[-iK-p(z)]>  must  be  computed  at  the  point  z'  =  z. 
Tatarski  has  used  variational  (or  functional)  derivatives  to  yield  a  result 
for  plane  waves [3],  but  it  appears  possible  to  obtain  similar  results  more 
simply: 

First  of  all,  we  note  from  Eq.(19)  that  p(z)  is  a  Gaussian  random 
variable  with  zero  mean.  Although  a  rigorous  proof  of  this  statement  is  not 
given,  it  can  be  verified  heuristically  by  applying  a  central-limit  theorem 
to  the  sum  of  sections  of  the  integral,  each  of  which  is  several  macroscales 
long. 

The  next  step  is  to  note  that 

<6e(K,z,)[p(z)]m>  -  m<6e(K,z,)p(z)>< [p(z) ]m-1>  (A.l) 

where  m  is  an  integer.  The  easiest  way  to  see  this,  is  to  write  out  the 

left-hand  side  as  a  product  of  m  integrals  dz1 - dzffi  after  substituting  Eq.(20) 

for  p(z).  The  correlation  is  then  essentially 

<6e(K,z')6e(K1,z;)) - 5e(Km,zm)>, 

a  product  of  m  +  1  factors.  Now,  whenever  z'  and  z^  (1  <  j  <  m)  are  both 
at  least  several  macroscales  L0  distant  from  each  and  everyone  of  the  m— 1 
remaining  z  coordinates  one  can  factorize  Eq.(A.2)  into 


(A. 2) 


<6e(K,z,)6e(Kj,Zj)><6e(K1,z1)' 


-/- 


(A.  3) 


where  the  slash  indicates  that  6e(Kj,Zj)  is  missing  in  the  second  correlation. 
Of  course,  it  is  also  true  that  we  could  have  separated  just  z'  and  then 
factorized  (A. 2)  into  <<Se(K,z')>  and  an  m-point  correlation,  but  <5e(K,z')>  =  0 
so  the  next  non— zero  factorization  to  look  at  is  Eq.(A.3)»  The  factorization, 
Eq. (A. 3) ,  is  not  possible  when  both  z'  and  z^  are  "near"  (with  respect  to 
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length  LQ)  to  one  or  more  of  the  remaining  m-1  points.  However,  that 

situation  restricts  the  dz, - dz  integration  to  an  effective  multiple- 

x  m  2 

integration  region  that  is  less  by  a  factor  of  order  (Lq/L)  than  in  the  case 
of  Eq. (A. 3) .  Of  course,  this  is  only  a  heuristic  description  of  a  theorem 
that  we  have  proved  more  rigorously  in  Sections  11  and  12  of  TRIII.  The  net 
result  is  that  we  may  replace  Eq.(A.2)  by 


P<Se  (K ,  z ' )  5e  (Kj  ,  z^ )  ><6e  (K1 ,  z^ - / - 5e(Km,  zm)  > 

It  is  easily  seen  that  Eq.(A.4)  yields  Eq. (A.l),  which  is  now  used  in 

the  following  sequence  of  steps: 

00 

<5e(K,z')e“i^'P(z)>  =  <5e(K,z')  [K- p(z)  ]m> 

m=0m- 


(A.  4) 


«Se(K,z')K*p(z)><[K*p(z)  ]m  1> 


(A.  5) 


=  <-i5e(K,z*)K"  p(z)><exp[-iK*p(z)  ]N 

2 

We  obtain  Eq.(21a)  from  this  last  step  by  noting  that  <exp£>  =  exp<£  /2> 
for  a  Gaussian  random  variable  £  with  zero  mean,  and  by  noting  that  x  and  y 
components  of  p(z)  are  equal  and  uncorrelated. 

^  ~y 

Now  we  derive  Eq. (21b)  by  substituting  Eq.(20)  into  <-i6e(K,z')K*p(z)> 


to  obtain 


Li 

-(8n2)"iy*dz1  M(L,z,z1) J c^KjK'K^e  (K,z,)5e(K1,z1)>, 

and  apply  Eq.(17)  of  TRV  [or  Eq.(5.2)  of  TRIII]  to  reduce  this  to 

L  L 

(e2K2/2)[(l-z/L)^z1z1$?(K,z'-z1)  +  ydz1(z-z1)$2(K,z,-z1)] 


(A.  6) 


(A.  7) 


where  $2(K,Az)  is  a  partial  spectrum,  namely  the  two-dimensional  Fourier 
transform  of  the  stationary  correlation  function 

<6e(p+Ap,z+Az)6e(p,z)>/<(6e) 2>  with  respect  to  Ap.  By  ignoring  some  boundary 

effects  of  0(L  / L)  we  can  reduce  Eq.(A.7)  and  obtain 
o 
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z'<z 


(A.  8) 


<-iSe(K,z')K*p(z)>  =  (e2K2/2)<J>(K)z’ (L-z)/L 

=  (e2K2/2)$(K)z(L-z')/L  ,  z’>z 

Note  that  the  two  results  in  Eq.(A.8)  become  equal  at  z'  =  z.  However, 

2  -* 

Eq.  (A. 8)  does  not  lead  to  the  right  answer  for  <6x  >  in  the  limit  p(z)  0. 

2 

The  limit  then  yields  a  value  for  <6y  >  that  exceeds  the  geometrical-optics 
approximation  by  exactly  a  factor  2.  Apparently,  there  is  a  factor  1/2 
missing  in  Eq.(A.8).  The  same  discrepancy  holds  in  the  plane-wave  case.  We 
have  included  an  extra  factor  1/2  in  Eq.(21b),  as  we  did  in  Section  3  of 
TRVI  for  plane  waves.  The  following  argument  is  suggested:  When  z'  >  z  it 

y  ^  -y  -y 

cannot  be,  physically,  that  p(z)  is  correlated  with  6e(K,z')  because  p(z)  is 

the  transverse  coordinate  of  a  ray  traveling  from  z^  =  0  to  z^  =  z,  and  the 

medium  at  z1  >  z  does  not  influence  the  medium  at  z^  <  z  beyond  several 

macroscales  L  .  It  is  true  that  the  boundary  condition  p(L)  =  0  forces  the 
°  ->- 
above  formal  result  for  z'  >  z  -  which  seems  to  imply  correlation  of  p(z)  with 

6e(K,z')  for  z'  >  z  -  but  the  ray  can  only  arrive  at  p(L)  =  0  because  the 

medium  beyond  z..  =  z  makes  the  ray  bend  into  p(L)  =0.  In  fact,  we  could 
-*•  -y  .  . 

physically  change  the  medium  beyond  =  z  and  it  appears  obvious  that  p(z) 

will  be  unchanged.  The  above  result,  Eq.(A.8)  for  z'  >  z,  appears  therefore 

to  apply  to  the  reciprocal  situation:  a  ray  traveling  from  z  =  L  to  z  -  0. 

This  situation  cannot  be  distinguished  in  this  formalism,  but  it  is  physically 

different.  Somehow,  the  information  about  the  direction  in  which  the  rays 

travel  has  been  lost.  If  a  time-dependent  formalism  were  pursued  more  carefully, 

we  would  obtain  a  Green's  function  in  ^(z),  hence  in  p(z),  that  would  be  zero 

for  "advanced"  times,  i.e.,  for  a  field  at  z1  =  z  due  to  a  signal  emanated 

from  z1  =  z'  when  z'  >  z,  and  that  will  rule  out  correlation  in  Eq.(A.8)  for 

z'  >  z  if  this  correlation  is  defined  more  carefully.  Note  that  reciprocity 

is  not  in  question:  the  result  for  z'  >  z  is  perfectly  symmetric  with  that 

for  z'  <  z  provided  the  ray  direction  is  considered.  The  suggested  form  of 

Eq. (A. 8)  for  rays  traveling  from  z^  =  0  to  z^  =  L  is 
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<-i<5e(K,z')K*p(z)>  «  (<2K2/2)<I'(K)z'  (L-z)/L  ,  z'  <  z  (A. 9) 

«0  »  z'  >  z 

and  a  reciprocal  version  "ith  zero  for  z '  <  z  when  rays  travel  in  the  other 
direction.  There  is  a  jump  discontinuity  at  z'  =  z  where  half  the  sum  of  the 
values  ior  z'  <  z  and  z'  >  z  in  Eq.(A.9)  should  be  taken.  That  finally  yields 

Eq. (21b). 


i 

I 
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